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Abstract. In a recent work Grappin et al. [1] have shown that low- frequency movements can be transmitted from one 
footpoint to the other along a magnetic loop, thus mimicking a friction effect of the corona on the photosphere, and invalidating 
the line-tying approximation. We consider here successively the effect of high frequencies and turbulent damping on the 
process. We use a very simple atmospheric model which allows to study analytically the laminar case, and to study the 
turbulent case both using simple phenomenological arguments and a more sophisticated turbulence model [2]. We find that, 
except when turbulent damping is such that all turbulence is damped during loop traversal, coupling still occurs between 
distant footpoints, and moreover the coronal field induced by photospheric movements saturates at finite values. 
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In large simulations the line-tying approximation is 
usually adopted for the coronal boundary. According to 
this approximation, strong density gradients cause the re- 
flection of any coronal disturbance. Hence, the photo- 
spheric velocity field can be prescribed since the reaction 
of the coronal dynamics is completely neglected. Apply- 
ing a velocity field uq at the base of a loop of length L 
and mean magnetic field Bo, keeping the other foot an- 
chored in the photosphere leads to an infinite accumu- 
lation of magnetic energy (driven by the velocity shear 
between the loop footpoints). To test whether the strati- 
fication really leads to line-tying, Grappin et al. [1] have 
studied the transmission of a signal from one footpoint 
to the other including a simple atmospheric model in 
a 1.5D MHD simulation. Footpoints are free to move, 
i.e. boundaries are transparent to waves, and a kick is 
given to left footpoint, injecting an Alfven waves from 
the left boundary. At the beginning, the magnetic en- 
ergy grows steadily according to the line-tied assump- 
tion, b ps B()Uot/L (transient acceleration of the left foot- 
point). On longer times T re ; = Ljvf 1 (V c ? h is the Alfven 
speed at the photosphere), the leakage of Alfven waves 
through the transition region (T.R.) accelerates (deceler- 
ates) the right (left) footpoint: the system relaxes to a 
state in which both footpoints have the same speed. 



SIMPLIFIED MODEL ATMOSPHERE 

We adopt a rough representation of the loop atmosphere 
modeling the transition region as a discontinuity in the 



Alfven speed, separating two uniform Alfven speed 
chromosphere and corona (hereafter subscript and/or su- 
perscript c, ph indicate quantities evaluated at the corona 
and at the photosphere, not distinguished from the chro- 
mosphere in the model). This allows to easily consider 
the effect of finite frequency, wave reflection, and wave 
damping while retaining the physical process responsi- 
ble for the relaxation. 

The small parameter e = v£ h / V a c quantifies the jump at 
the T.R. and hence reflection. Continuity of magnetic 
and velocity field fluctuation, 8b, 8u, at the two T.R.s 
yields four jump conditions for the Elsasser variables 



= <5m=f 8b/y/4np: 



-u L , u\ — U\ — {u^ — u L )/e, 
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u\-u R =u\ jz, (m 3 =0), (2) 



here labeled uf LR y Superscripts ± indicate the di- 
rection of propagation (rightward and leftward respec- 
tively), while 1,L,R,3 specifies the position with re- 
spect to the T.R. where they are evaluated: left chromo- 
sphere, left corona, right corona, right chromosphere re- 
spectively. We further use the notation ph, 1 and ph, 3 for 
evaluations at the left and right footpoints (see Fig. 1). 
We chose solutions of the form u °< exp[— icot — t/x\ 
in which co is the wave frequency and T is a damping 
timescale. A wave of amplitude u\ j = u® is injected 

at the left footpoint so that u\ = w^expf— i(Ot p h\. Time 
is counted in coronal crossing time t = nt c = nL c /V£, n 
indicating the number of reflections, so that wave prop- 
agation and damping in the corona (we exclude damp- 
ing in the chromosphere for reasons that will be clear 
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FIGURE 1. Modeled loop atmosphere (top) and jump con- 
ditions (red circles) at the two T.R.s (bottom) for the rightward 
and leftward propagating waves (u respectively). 



after) are accounted by the factors | = exp[— i(Ot c ) and 
J3 = exp[-* c /r]: 



it 
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= ZPut\„, m l|„+1 =ZP u R \ n - 



(3) 



Combining Eqs. 1-3 during a cycle of reflection and 
propagation in the corona yields a recurrence formula, 

fl=(l -£)/(!+£), 



u'lln = (1 +a)u^\n + a 2 li 2 ^ 2 u^\ n -2 



(4) 



(and a similar expression for m^) that can be used to 
evaluate for each variable uf: a) its asymptotic value, 
setting n = n — 2; b) its temporal evolution, expliciting 
«U-2 = /("i1n-4, «|n-4) and so on. 
Here we consider only two special cases: a) Finite fre- 
quency fluctuations ^ 1) with no damping; b) Zero 
frequency fluctuation with damping (j3 ^ 1). 



NO DAMPING, co^O 

With j3 = 1 , the asymptotic solutions for the reflected and 
transmitted wave can be written as 



(l-« 2 )^ 



•(l-a 2 ! 2 )" 1 ' "' 3 (1-fl 2 ^ 2 ) n 
while in the corona one finds 



(l-a 2 <«) 

(l+«)g S 
(1-^ 



(5) 

(6) 
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For injected waves of vanishing frequency and very low 
coronal density (small e) one finds immediately from 
Eqs. 5-7 that uf w 1/(1 — a)uf w 1/euf ^> and = 
(1 — a)m^ « 2euf. Hence full transmission («+ « M3") 
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FIGURE 2. Asymptotic shear between footpoints Am (bot- 
tom) and coronal magnetic field (top) as a function of fre- 
quency. 



implies very large coronal amplitudes uf s» 1/e (larger 
than the WKB prediction u WKB ml/ ^e) and also uf w 
— m^. The latter finally implies db/^fAnp ^s> 5«, that is 
magnetic energy much larger than kinetic energy in the 
corona. The asymptotic solution is actually a uniform 8u 
and 8b along the loop, while line-tying would have given 
a linear profile for 8u (from u® to 0) and 8b that goes to 
°° (never ending shear) 1 . 

The asymptotic velocity shear between the right and 
left footpoints Au = \8u p h\\ — |5m ; ,/, 3| and the coronal 
magnetic field 8b c for finite frequency fluctuations are 
shown in Fig. 2 as obtained from Eqs. 5-7 assuming 
a 10 Mm loop, L ph = 2 Mm, Vf 1 ' = 0.7 km/s, and 
£ = 1 /70 <C 1. 8b c ranges from the photospheric value, 
8b „h = 15 Gauss, at zero frequency to zero at high 
frequencies. Because of the almost perfect reflection at 
the two T.R.s the coronal part of the loop behaves as 
a resonant cavity (compare the poster [3] in the same 
session and [4]), producing peaks at the fundamental 
frequency V/j /L c and its harmonics. On the other hand 
the asymptotic shear (Fig. 2, bottom panel) has well 
defined peaks well below such frequency that are caused 
by the phase shift between the injected and reflected 
waves at the left footpoint, i.e. a dependence of the form 

Au = g(03t p h), t p h — Lph/Vcf' 1 being the chromospheric 
crossing time. Figure 3 shows the temporal evolution of 
Au for different frequencies (black lines) for the same 
10 Mm loop. The relaxation is independent of e (cf. the 
plots for co = Hz obtained with e = 1 /70 and e = 0.4), 
which allows to derive a compact formula for Au(t) by 
expliciting the dependence m„_2 = /(m„-4) = f{ u n-(>)--- 
and taking the limit e — > 0: 



= ufa[l-{l-a 2 )^ 2 J^(t)] (8) 
= uf(l-a 2 )E,.Jf(t) (9) 



the line-tying can be recovered taking the limit e — > before n 
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FIGURE 3. Velocity difference between footpoints Am as a 
function of time for different frequencies, for CO = and e = 0.4 
(blue-dashed), and for the weakly nonlinear case (red-dotted). 



with; jT(f) 



1 — exp [— 2t I Xrel - icot] 
1 + i(DT re i 



(10) 



and irei — L c /Va h being the relaxation timescale analo- 
gous 2 to that found by Grappin et al. [1], T re / = L/V£ h . 
The peak at t w 2f p /, corresponds to the initial accelera- 
tion of the left footpoint exerted by the first reflection at 
the T.R. (transient line-tying). The following relaxation 
is due to the deceleration of the left footpoint and the ac- 
celeration of the right footpoint. Oscillations in Au(t) for 
finite frequencies are caused by the phase shift between 
the injected and reflected waves at the left photosphere 
(the icot term). 



TURBULENT DAMPING, (0 = 

We are now interested in understanding whether the 
above relaxation (for zero frequency) takes place in pres- 
ence of a damping, that we chose to be caused by incom- 
pressible turbulence triggered by counter-propagating 
waves. We perform simulations in which the nonlin- 
ear MHD dynamics in planes perpendicular to the mean 
magnetic field is replaced by two-dimensional shell mod- 
els for 8u and 8b. The shell models of different planes 
are coupled by Alfven waves propagating along the mean 
magnetic field [5, 2]. In each layer the damping factor 



is given (dimensionally) by /3 = exp[ 



ph, c 



t p h. 



ph, 



with 



/ 



i i / U„, aiiu / 
taken to coincide with the loop width. We expect that 
the amplitude in the corona saturates to value smaller 
than u^r oe 1/g when subject to a damping, say the WKB 
value u c «* i/y/e. One gets t c /r^ = (u° iy /e/Va ph )(L c /l ± ) 
that justifies neglecting the damping in the chromo- 
sphere. Based on this dimensional analysis, we set u\ = 



and / 1 is the outer scale of turbulence 



0.2 km/s, L c = 6 Mm e = 0.4, L c /l ± = 1, u\jvf =2/1 
in the simulation to obtain a weak turbulence limit 3 . The 
velocity shear between footpoints relaxes in a few char- 
acteristic timescales, although a small shear survives, 
oscillating with time (red-dotted curve in Fig. 3). For 
smaller £, relaxation occurs earlier and to larger Am 
(cf. symbols in Fig. 4, A), with oscillation on similar 
timescale (not shown). The nature of the oscillations is 
not completely clear, although it is related to the sharp 
reflection at the T.R. and is expected to be reduced when 
the T.R. is properly accounted for. 
The asymptotic shear can be understood in our simple 
analytical picture adopting a phenomenological expres- 
sion for the turbulent damping. With j3 < 1 and | = 1 in 
the recurrence formula, and in the limit of small £ one 
obtains for the wave amplitude in the corona: 



1 1 -exp[-2f/T re /-f/T„/] 



2e 



1 + T re/ /2v 



(11) 



At the beginning, when few reflections have occurred, 
Z„l ^> z re i and the system evolves as in the linear case, 
with the wave amplitude in the corona growing toward 
the asymptotic non-WKB value u[ =u { \/e. When T„/ < 
Z re l the relaxation process stops and the level of fluctua- 
tion saturates to the value: 



1 



1 



2g 1 + XreljlXrd 



(12) 



Assuming that u) 



(or equivalently 8u <C 



8b c / ' y/Aizpc), and choosing a Kolmogorov type phe- 
nomenology for the nonlinear timescale [6], V = l±/ u^, 
yields the asymptotic level of the coronal magnetic field 



8K 



\JAnpc 



T pf> 




(13) 



Given the value of m+ w m^, the velocity shear at foot- 
points can be obtained from the asymptotic expression 
or the time evolution formula for Am. Here we prefer to 
exploit the stationary state of the system in the induction 
equation, which in a dimensional form leads to (other 
choices gives similar results): 



,. Am bee 

v c 

a L c 



0, 



from which one finally gets, 



Am = 



VI 



P h 



(14) 



(15) 



2 In their model atmosphere, Grappin et al. [1] have non-vanishing 3 Keeping fixed the ratio u^/Va 1 ', a different choice of the nonlinear 
Alfven speed gradients at the photospheres, corresponding to our T.R.. timescale, as done below, would give a slightly higher e. 
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FIGURE 4. Velocity shear Au (A) and asymptotic coro- 
nal magnetic field Sb c (B) as a function of e from eq. 12 
and eq. 14 solved with different nonlinear timescales: a Kol- 
mogorov (dashed) and an IK (solid) type. Symbols are results 
from the shell simulations. 



With a different choice of the nonlinear timescale, as in 
the IK phenomenology [7, 8], T„/ = I jV£ / \u^) 2 , one 
obtains a cubic expression for m+, which has to be 
solved numerically, but the induction equation can still 
be used to get the velocity shear. Note that for small e, 
Eq. 13 gives a dependence consistent with the WKB es- 
timate u^r °c £~ 1//2 , while the IK phenomenology gives 
uf oc g~ 2 / 3 , that is a coronal magnetic field intermedi- 
ate between the WKB and linear solutions. In Fig. 4 
symbols represent 8b c and Au computed at the coro- 
nal base obtained with the simulations (and averaging 
over 4-5 oscillations), varying e while keeping fixed the 
other parameters (loop length and width, Alfven speed 
at the photosphere, injected wave amplitude). The solid 
and dashed lines are the analytical estimates obtained 
with a Kolmogorov and IK phenomenology. A relatively 
good agreement is found between the simulations at high 
Alfven speed contrast (small e) and the IK phenomenol- 
ogy. For the latter a fiarly accurate approximation is 
given by (cf. Fig. 4, B for e < 1(T 3 ): 



This situation continues to hold basically for CO <C ft>* ~ 
1 / % re \. For CO « ft)*, the asymptotic shear is finite (partial 
relaxation) and 8b c < 8b „/,. For CO ^S> ft)*, our model cor- 
rectly reproduces the behavior of the coronal magnetic 
field that tends to and shows resonances at well defined 
frequencies. On the countrary oscillations in the veloc- 
ity shear, < Au = Mj, are exaggerated in our simplified 
model but when the T.R. is represented as a continuous 
layer we recover the correct behavior Am — > 0. The sys- 
tem is neither line-tied (8u = 0) nor relaxed (8b c = 0). 
When weak (turbulent) damping is taken into account, 
both numerical simulations and analytical calculations 
show that the relaxation process still occurs. In the strong 
turbulence limit, the loop relaxation is stopped as soon 
as the coronal nonlinear time becomes smaller than the 
coronal propagation time (if, < t c ). The asymptotic mag- 
netic field in the corona (in velocity units) has a scaling 
that varies between e -1 / 2 and £~ 2 / 3 , depending on the 
phenomenology (Kolmogorov and IK respectively, the 
latter showing a better agreement with the simulations), 
which is intermediate between the linear (1/e) and WKB 
(l/Ve) scalings. The relaxation still occurs when the 
ratio TS/T re / > 1. According to the definition of non- 
linear timescale in the Kolmogorov or IK phenomenol- 
ogy, this can be rewritten as x L nl jz re i = (erff /r re i) 1 / 2 

and xf,/T re / = (er^/Tre/) 1 / 3 respectively. For a given 
loop width and initial footpoint shear, and assuming that 
longer loop has higher density contrast (smaller e) the 
asymptotic state is controlled by loop length. 
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{8b° ph — 1 /2u\yjAnp p h ~ 15 Gauss). For small reflection 
(e > 0.02), m+ u~, 8b c /y/4np c no uf, so that both phe- 
nomenologies fail to reproduce the simulation results. 



DISCUSSION 

Adopting a simplified model based on Alfven wave prop- 
agation and reflection we recover the essential features 
found in Grappin et al. [ 1 ] . A loop driven at one footpoint 
by a "zero-frequency" photospheric velocity relaxes to a 
uniform transverse velocity Au = and magnetic field 
8b c = 8b p h on a characteristic timescale T re / = L c /V£ h . 
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